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ABSTRACT 

In order to understand the collapse dynamics of observed low-mass starless 
cores, we revise the conventional stability condition of hydrostatic Bonnor-Ebert 
spheres to take internal motions into account. Because observed starless cores 
resemble Bonnor-Ebert density structures, the stability and dynamics of the star- 
less cores are frequently analyzed by comparing to the conventional stability 
condition of a hydrostatic Bonnor-Ebert sphere. However, starless cores are not 
hydrostatic but have observed internal motions. In this study, we take gaseous 
spheres with a homologous internal velocity field and derive stability conditions 
of the spheres utilizing a virial analysis. We propose two limiting models of spon- 
taneous gravitational collapse: the collapse of critical Bonnor-Ebert spheres and 
uniform density spheres. The collapse of these two limiting models are intended 
to provide the lower and the upper limits, respectively, of the infall speeds for 
a given density structure. The results of our study suggest that the stability 
condition sensitively depends on internal motions. A homologous inward motion 
with a transonic speed can reduce the critical size compared to the static Bonnor- 
Ebert sphere by more than a factor of two. As an application of the two limiting 
models of spontaneous gravitational collapse, we compare the density structures 
and infall speeds of the observed starless cores L63, L1544, L1689B, and L694-2 
to the two limiting models. L1689B and L694-2 seem to have been perturbed to 
result in faster infall motions than for spontaneous gravitational collapse. 

Subject headings: Star: formation — ISM: kinematics and dynamics — ISM: L694-2, 
L1689B, L1544, L63 
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1. INTRODUCTION 

Low-mass star formation is a fundamental process in astrophysics because the majority 
of stars in our Milky Way are low- mass stars. Recent observations at (sub) millimeter and 
IR wavelengths show that starless cores are the initial stage of low-mass star formation 
(Myers & Benson 1983; Benson & Myers 1989; Ward- Thompson et al. 1994). Isolated 
low-mass starless cores represent the simplest environment of star formation compared to 
high-mass star forming region (e.g. Orion, Cygnus, Carina, etc.) because they are in less 
crowded regions with less contamination from strong radiation fields of neighboring stars 
and from feedbacks of stellar outfiows. A starless core usually forms a single star or a few 
stars at most (e.g. Lada 2006). Starless cores are therefore ideal objects to research the 
basic physics of star formation. 

Star formation involves a series of complicated dynamical, chemical, and radiative 
processes. In one limit of theoretical studies, starless cores are formed through gravitational 
fragmentation and undergo hydrostatic collapse when they become gravitationally critical 
through loosing their thermal or magnetic pressure support (e.g. Mouschovias & Ciolek 
1999). In this limit, gravitational contraction is the dominant mechanism during collapse. 
Another limit of theoretical studies, which is still developing, proposes that turbulence 
is comparably important to self-gravity in forming stars (e.g. Mac Low & Klessen 2004; 
Kudoh & Basu 2008; Nakamura & Li 2008; Gong & Ostriker 2009; Gong & Ostriker 2011). 
In this limit, colliding fiows in supersonic turbulence create or disrupt starless cores, and 
some fraction of starless cores made from colliding fiows may undergo collapse. The main 
differences between this limit and the former one are that the speed of colliding fiows 
may be much faster than the speed of gravitational contraction and that the ram pressure 
of colliding fiows may considerably shorten the time to form stars. Thus, studies of the 
internal velocity fields of starless cores are necessary to determine which limit is dominant. 
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The internal velocity fields of starless cores are probed using dense gas tracers such as 
HCO+, HCN, CS, etc. (e.g Lee et al. 1999, 2004; Gregerson et al. 2000; Evans et al. 2005; 
Sohn et al. 2007; Stabler & Yen 2009). Transitions of these tracers are usually optically 
thick lines which have an asymmetric profile compared to optically thin symmetric lines 
when there are internal motions. Surveys of the skewness of CS 2—1, CS 3—2, and HCN 
1-0 lines were done by Lee et al. (1999, 2004), and Sohn et al. (2007). Their studies 
show that 25%, 27%, and 20% of their samples, respectively, have blue-skewed line profiles 
indicative of infall. In addition, studies on the internal velocity fields of L694-2 and LI 197 
with radiative transfer models and dynamical simulations (Lee et al. 2007; Seo et al. 2011b) 
suggest that both starless cores have supersonic infall motions substantially faster and 
occurring in a narrower layer than predicted by quasi-static gravitational collapse. Such fast 
infall motions may appear in a scenario of colliding shocks where the turbulent surrounding 
may enhance infall speed in a narrow layer (Gomez et al. 2007) and considerably affect 
stability of the starless cores because the ram pressure Pram ~ pv"^ becomes comparable to 
the surface pressure Psurf ~ pc^. These studies demonstrate that internal motions of some 
starless cores may not be owing completely to gravitational contraction, but may also be 
under considerable influences from their surrounding environments. Their internal motions 
must be taken into account when testing their stability and understanding their internal 
dynamics. 

Dynamical evolution of starless cores in the limit of negligible external dynamics 
(gravity dominant limit) are relatively well studied utilizing similarity solutions and 
numerical simulations (e.g. Larson 1969; Pension 1969; Shu 1977; Hunter 1977; Foster & 
Chevalier 1993; Aikawa et al 2005; Seo et al. 2011b). The fiducial model in this limit is the 
collapse of a self-gravitating, isothermal, hydrostatic gas sphere bound by external pressure 
called the Bonnor-Ebert sphere (Ebert 1955; Bonnor 1956; Ebert 1957; hereafter BE). 
This model has been widely used because it is a physically well established model with its 
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column density profile similar to observed column densities in near-IR absorptions and dust 
continuum emissions (e.g., Alves et al. 2001; Evans et al. 2001; Tafalla et al. 2002; Harvey 
et al. 2003a, 2003b; Kandori et al. 2005). In this model, the critical BE sphere collapses 
while maintaining a BE-like density profile (Foster & Chevalier 1993; Aikawa et al 2005; 
Seo et al. 2011b). However, the initial BE sphere is a hydrostatic sphere neglecting any 
internal motion, whereas real starless cores are likely contracting, expanding, or oscillating 
(e.g. Lee et al 1999, 2004; Keto et al. 2006). The dynamics of starless cores in the limit 
of external dynamical influences (turbulence regulating limit) are also studied through 
simulations (e.g. Kudoh & Basu 2008; Nakamura & Li 2008), but do not have a simple 
fiducial model because of the chaotic nature of turbulence. Since only the first limit has a 
fiducial model, it is hard to compare observed starless cores to the simulations and interpret 
their dynamics with respect to the two limits. A dynamics model is required that provides 
a boundary between the two limits. 

The main purpose of this work is to take an internal velocity field into account in the 
determination of the dynamical stability condition for starless cores. Our stability condition 
serves as a more realistic stability criterion for starless cores with internal motions. We 
also suggest a range of infall velocities at a given density structure that may classify 
observed starless cores with respect to the two dynamical evolution limits. The comparison 
of observed infall speeds of starless cores to this range will tell us whether the collapse of 
starless cores mostly depend on its internal properties or evolve together with its turbulent 
surroundings. 

In this paper, first, we derive a new dynamical stability condition of the BE spheres 
with various dimensionless size and homologous infall motions. Hunter (1979) studies 
stability of a uniform field with homologous infall motions using the virial treatment. We 
apply Hunter's analysis to the BE spheres and derive critical sizes of BE spheres for a given 
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homologous infall motion. We define ^max as the dimensionless size 

Smax — 5 \^ ) 

Cs 

where t is time, and G is the gravitational constant, Rs and pc are the outer radius and the 
central density of the core. We calculate dynamical domains in a two dimensional ^max "^s. 
peak infall velocity space. The critical size ^crit for a given infall speed will separate the 
stable and unstable regimes in the ^max vs. peak infall velocity space. 

Second, we search a range of infall velocities for isothermal spheres without any initial 
internal motion and any external perturbation. This is to find dynamical domain of the 
pure gravitational collapse in the ^max vs. peak infall velocity space. Without any initial 
internal motion and external perturbation, the dynamics of the starless core depend only 
on the density distribution and there are two limits of the density distributions for starless 
cores: the BE sphere and a uniform sphere. The collapse of the critical BE sphere gives the 
slowest infall with respect to the density concentration (Seo et al. 2011b) and provides the 
lower limit of the infall speed in spontaneous gravitational collapse (or pure gravitational 
collapse). A starless core that collapses slower than this model is likely to be perturbed 
by external perturbations or is not isothermal. The uniform sphere is a starless core with 
the least pressure gradient against the self-gravity for its initial conditions. Inevitably, the 
sphere develops a faster infall motion than the collapse of the critical BE sphere at a given 
density structure (Seo et al. 2011b). Since the uniform density field has the shallowest 
profile that a starless core can have, the collapse of the critical uniform sphere has the 
fastest infall speed with respect to a given density concentration (in the absence of external 
perturbations). Any starless cores collapsing faster than the collapse of the critical uniform 
sphere is likely to have initial converging flows, or be perturbed to collapse faster than any 
pure gravitational collapse by an external source. 

As an application of our study, four starless cores (L1544, L1689B, L63, and L694-2) 



-7- 

are compared with the new stabihty condition and the domains of two dynamical hmits 
(gravity dominant and turbulence regulating limits). The four starless cores are the best 
four infall candidates with studies of both density structures and internal motions in dust 
continuum or extinction and molecular tracers of infall (e.g. HCN, CS, and HCO'^). In 
order to normalize the physical quantities of the four cores into dimensionless quantities 
of this study, we need to measure gas temperatures of the starless cores. We observed the 
four starless cores in NH3 (1,1) and (2,2) inversion lines with the Robert C Byrd Green 
Bank Telescope and deduced their gas kinetic temperature. Their collapse dynamics are 
interpreted with the results of our dynamical calculations. 

The layout of this paper is as follows: §2 briefly introduces our observation of NH3 lines 
measuring temperatures of the four starless cores. In §3 we introduce the result of Hunter 
(1979) and elaborate our application of his method to an arbitrary density distribution. In 
§4 we outline the two limiting models of the spontaneous gravitational collapse, and discuss 
the physical meaning of the regimes divided by the two models in the ^max vs. peak infall 
speed space. In §5 we discuss limitations of our analysis and results. We also compare 
observations of the four starless cores with our stability conditions and discuss their internal 
dynamics. Finally, in §6, we will summarize our results and discussion. 

2. GBT OBSERVATIONS 

Observations of the NH3 (1,1) and (2,2) inversion transitions were performed over three 
observing shifts in 2006 (GBT06A68) with the Robert C Byrd Green Bank Telescope. 
Each observations consisted of a 4 minute frequency switched observation with a 4.11 MHz 
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throw at 4 Hz centered on the 850 /xm continuum peak position (see Table 1). The GET 
spectrometer was set up with 50 MHz bandwidth and 6.1 kHz spectral resolution (0.08 
km/s resolution). The spectra were folded and baselined using standard GBTIDL routines. 
The main beam efficiency was determined from observations of quasars (3C48 and 3C286) 
and planets to be r]mb = 0.74 ± 0.05. In Section 5 we analyze the NH3 observations for the 
four starless cores with measured infall motions (L1544, L1689B, L63, and L694-2). 

3. GRAVITATIONALLY STABLE OR UNSTABLE: VIRIAL TREATMENT 

To facilitate the virial analysis for a starless core, we assume the core is an isothermal 
gas sphere bounded by external pressure. We also assume that all internal motions 
are spherically symmetric. The virial theorem for an isothermal medium with spherical 
symmetry is given by 

1 M rRs rRs rRs 

-—^ / p{r)r'^ ■ Airr'^dr = / A'nr'^p{r)r'^dr + 3 / A7rr'^p{r)c1dr 

2 dt Jq Jq Jq 

-/ Anr^pir)^^ ^-^^ — rfr - 47ri?^c^p(i?,), (2) 

Jo '^ 

where Rg is the radius of the core, and Cg is the sound speed. The above equation may be 
normalized with the following normalizing parameters: The sound speed at the core is 

/ T \^/^ 
c,,e = 0.188 km s-i (tok) + ''°*' ^^^ 

where T means the temperature of the core, and a^t is the rms speed of non-thermal 
components. The normalizing parameter for length a is defined as 

" " VSUK = °°** "' ( o,188k,ns-0 (^^j (— ) • <^' 

where Uc is the number density of molecules at the core center, and p is the mean molecular 
weight. For the normalization of time, we take the sound crossing time to) which is the time 
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taken by the sound wave to travel a unit length, 

1 , /lO^cm-^V^' /2.33\'/' 
to = -jr^= = 0.23 X 10^ yr . 5 

The normalized form of equation ([2]) is 

-2s(es)4', (6) 

where s = p/p^t = 0), ^ = r/a, ^s = Rs/a, M(0 = Jq ^^s{^)d^ and r = t/to- 



3.1. A Uniform Density Core with a Homologous Motion 

The dynamical stability condition of a uniform density field with homologous infall 
motions is derived by Hunter (1979). He calculated how fast the converging motion should 
be to make a stable core into a critical configuration after contraction. In this paper, we 
show his results and briefly discuss its meaning in ^^ax vs. peak infall speed space. For 
detailed derivation, please refer to his paper. 

The relation between the initial speed of homologous motion and the ratio of initial 
and final size of the core (equation (4) in Hunter's paper) in terms of our dimensional 
variables is as follow: 

/a 1 \ / 1 \ " 

(7) 

where Ui is the initial peak velocity of homologous motion divided by the sound speed. R 
is the compression ratio defined as ^si/^sf where ^si is the core radius at the initial moment 
and ^sf is the core radius at the final moment. When there is no external pressure, equation 
([7]) becomes 

U? = 10 ( lnR + — -1 
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The initial velocity f/; are plotted in black lines as functions of the compression ratio 
R in Figure [TJ The solution oi U^ < is not a true solution because the velocity becomes 
imaginary. The solution of f/^ > and R > 1 is for a collapsing core, while the one of 
R < 1 is for a expanding core. The solutions for the cores bound with the external pressure 
of Pext = c^'S(.^si) and Pext = are marked with the solid and dashed lines, respectively. 
The solution of Pext = cisl^si) has imaginary values from R = 1 to 2.25. This is because 
we assume that the density structure is kept as a uniform field during the contraction, 
which is unrealistic. Because of this assumption, the work done by the external pressure 
during the contraction cannot be correctly estimated. On the other hand, the solution 
without the external pressure has no imaginary values since the stability depends only on 
the gravitational potential energy and the kinetic energy. 

Since core mass is fixed and we assume the uniform density structure is maintained, 
the radius and the central density of the core after contraction are Ps/R and pc{t = 
t{) = Il^pc{t = 0), where t{ is the time epoch at the completion of contraction by inward 
motions. So, the dimensionless size ^max becomes ^max(i = U) = ^maxit = 0)R^''^. The final 
dimensionless size ^max(^ = tf) of the core should be the same as the critical size of the 
static uniform density core. The critical size of the pressure-bound, static uniform density 
core is ^max =2.25, which is numerically estimated using a Godunov-type hydrodynamics 
code within an error of 5% (Seo et al. 2011a, 2011b). For a given inward velocity Ui and 
the final critical size .^max(^ = ^f), the compression ratio R is determined, and we can derive 
the corresponding ^max(i = 0) using equations (171) and (E]). The stability diagram for a 
uniform density cores with inward homologous motions is presented in Figure [21 The black 
solid line represents the core confined by the external pressure, and the black dashed line 
represents core with no external pressure. The regime at the left side of the black dashed 
line is where a core becomes gravitationally stable, while the one at the right side of the 
black solid line is the unstable domain for the critical uniform density core. The stability in 
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the regime between the two hnes depends on the external pressure. 

3.2. The Bonnor-Ebert Sphere with a Homologous Motion 

Observed starless cores are not a uniform density cores, but are reported to have BE-like 
density profiles. If the density field is not constant, but a function of the radial distance, 
the moment of inertia / and the gravitational energy W take forms of Co{^s)M{^s)^s ^''^^ 
Ci{C,s)M'{C^s)/(,s, respectively, where Co(^s) and Ci{^s) are functions of the radial distance. 
If the Cq and Ci do not sensitively depend on the radial distance and time during collapse, 
we may assume them as constants. The Cq value decreases from 0.2 to 0.14 for .^max from 
to 7 and becomes almost constant for ,^max >10. Since the total variation of Cq is about 
30% for smaller cores of ^max < 6.5 while the observational uncertainty in ^max is at least 
a factor of two (Kandori et al. 2005), then assuming Cq as a constant is a reasonable 
assumption unless the core is considerably disturbed and deviates from a BE-like structure 
(quasi-equilibrium contraction). Likewise, the Ci value varies from 0.6 to 1 (40%), and thus 
we may also assume it is approximately a constant. In this study, we take Co and Ci as 
Co = C'o(Csi) and Ci = Ci(^si), respectively. 



The virial theorem for the BE sphere is given by 

,CoM(6)es = 2i'CoM{Q + 2M{Q - 2Ci^-^-^ - 2eAQ- (9) 



—CoMi^es = 26 CoM{Q + 2M{Q - 2Ci— ^ - 2e,s,^sj 



We assume that the core does not get any additional mass during contraction so the core 
mass Mo is a constant Mo = M(6i). Assuming that the inward motion is a homologous 
motion, we have 

^^ CoCs CoC' Co Mo' ^ ^ 

We integrate the above equation from ^f to ^i and assume that the external pressure is 
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constant s(^si) = s{C,s{). Then, 



r2 I 2 i„^^sM I o^l A~:r I ^ ^\ ^ ^i^sijf^s ^3 



U-f = Ut + ^ In(^) + 2-iM. ^ - r - ^^(^si - e^f), (11) 

(--0 ^sf <-^0 \?si ^sf/ JL^O Mo 

where Uf is the velocity of the homologous motion at the final moment. Since the total 
mass of the core is fixed, we may write the core mass Mo = C2(^si)^si- The equation flTTl) 
becomes 



2 , esi , ^c^iMo /Csf ,\ 2 ,^ , / e; 



We assume that the kinetic energy of the homologous motion is expended at the final 
moment (f/f = 0) and the core becomes relaxed (^sf = 0). Equation ( 1T0|) at the final moment 
is 

1 gfcOeff _CiMo 
Putting the above result into the equation (IT^ . we get 
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U? 
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(14) 



When there is no external pressure, the above equation becomes 



In Figure [1] the square of initial velocity of equations f lT^ and flTSj) are plotted in blue 
lines as a function of the compression ratio R. The solid blue line represents the solution 
of Pext = c^s(^si), and the dashed blue fine represents the solution without the external 
pressure. Only the solutions at t/j^ > and R > 1 correspond to collapsing cores. In 
this figure, the Co is fixed to be Co(^max = 6.5). The solution of Pext = c^'5(^si) also has 
imaginary velocity at the range of R from 1 to 1.5 because the assumption of keeping a 
BE-like structure during the contraction results in incorrectly estimating the work done by 
the external pressure. 
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The critical size of BE cores with inward homologous motions can be calculated as 
follows: The final radius of the core after contraction by the inward motion, Rgf, is given by 

-Rsf = tt^sf = ttf^max.f, (16) 

where at and .^max,f are the newly defined normalizing parameter and dimensionless size of 
the core, respectively, after the contraction. Since a ^^ pc , the equation (ITBI) may be 
reduced to 



1 



where the pd and pd are the central densities at the initial moment and after the contraction, 
respectively. The relation between the central densities is given by the conservation of the 
total mass of the core, 
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where ip = —Ins. Reducing the above equation and combining with the equation fITTj) . we 
have 



R 



Smax,t j^f 






(19) 
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To be a critical core, the final size of the core ^max,f should be equal to the critical size of 
the static BE sphere, ^max,f = 6.5. If the initial size is given, the compression ratio R can 
be estimated for a critical core using the above equation. The corresponding peak velocity 
of the homologous motion can be calculated from the equation ( 1T5|) . Thus, we may relate 
an internal motion that makes a core to be critical with an initial size of a BE sphere. 

The critical sizes of BE cores with homologous internal motions are plotted in Figure 
m The blue solid and dashed lines are the critical sizes of the BE cores with the external 
pressure of Pcxt = ^^^(^si) and without external pressure, respectively. Because of the 
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imaginary velocity in the solution of Pext = c^,s(^si), the ^max value of the blue solid line 
starts from much smaller values than the critical size of the static BE core. The choice of 
Cq and Ci changes the starting points of lines for pressure bound cores, but not by more 
than factor of two. The critical sizes of smaller BE cores converge to those of uniform 
density cores because smaller BE cores have flatter density profiles. The critical size is very 
sensitive to internal motions. A homologous inward motion with a transonic speed can 
reduce the critical sizes of uniform density cores and BE spheres by more than half. 

4. SPONTANEOUS COLLAPSE OF STARLESS CORES AND LOWER 
AND UPPER LIMITS OF INFALL SPEEDS 

The collapse of the critical BE spheres and uniform density spheres provide the lower 
and upper limits, respectively, of infall speeds in the spontaneous gravitational collapse 
of isothermal cores. Figures 3 and 4 of Seo et al. (2011b) show density and velocity 
distributions of the collapse of the critical BE and uniform density spheres, respectively. 
Their density structures resemble each other at the later stage of evolution, but infall 
motions of the collapse of the critical uniform density sphere is always faster than the 
collapse of the critical BE sphere. 

From the numerical simulations of the collapse of the critical BE and uniform density 
spheres, we measure the peak infall velocities of the two models as a function of ^max- The 
results are plotted in Figure [31 The red sohd line represents the collapse of the critical 
BE sphere, and the red dashed line represents the collapse of the critical uniform density 
sphere. The stability lines are also plotted in the same figure with the same colors and 
line styles used in Figure [2J The collapse of the critical BE sphere shows a steady and 
slow growth of infall speed along with ^max- On the other hand, the collapse of the critical 
uniform density sphere results in a sudden increase of infall speed without a considerable 
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density concentration because there is no pressure gradient at the initial moment. Owing 
to the sudden infall motion at the early evolution stage, the core bounces momentarily 
and generates a strong accretion shock, which travels outward and regulates the growth 
rate of infall motion. The accretion shock is shown as a kink or cusp in the velocity fields 
in Figures 3 and 4 of Seo et al. (2011b). The kinks in Figure [3] at ,^max = 8 and 30 are 
when the accretion shock passes through the peak infall layer and when it escapes the core, 
respectively. 

Since the collapse of the critical BE sphere provides the lower infall speed limit of 
the spontaneous gravitational collapse, the regime above the red solid line implies that 
the gravitational collapse is hindered by an external perturbation. The regime below the 
collapse of the critical uniform density sphere is for gravitational collapse enhanced by an 
external perturbation. For a starless core that evolves between the two lines, it is hard 
to tell which factors, either initial density distribution or external perturbation or both, 
determines its dynamical evolution. Both factors may be equally responsible. 

5. DISCUSSIONS 

5.1. Limitations of Our Analysis 

The stability of isothermal gas spheres with homologous motions is derived utilizing the 
virial theorem. We impose homologous internal motions because they facilitate an analytic 
analysis. However, there is no observational evidence that starless cores have homologous 
internal motions when they are collapsing or perturbed. Moreover, a starless core does not 
maintain a homologous velocity field during its collapse even when a homologous velocity 
field is given as an initial perturbation. Nevertheless, our virial analysis provides a rough 
estimate of how much the conventional stability condition of BE sphere varies when there 
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are inward velocity fields. Results in this paper suggest that the stability condition for 
starless cores is very sensitive to internal motions. A transonic homologous inward motion 
can reduce the critical size of the static uniform density and BE cores by about 67% and 
54%, respectively. 

Spherically symmetric starless cores are assumed in this study. For a spheroidal core 
in a hydrostatic state, its critical size is the same as the critical size of the BE sphere as 
long as its projected area is the same as that of the BE sphere (Lombardi & Bertin 2001). 
Since the virial theorem is essentially an energy analysis, a stability analysis for spheroidal 
cores would be the same as this study except for the values of Cq and Ci. Infall motions 
of the non-spherical cores were studied by Myers (2005) whom showed that the peak infall 
speeds are about 50% faster for a cylindrical core and about twice faster for an infinitely 
extended slab. Observed starless cores usually have <2:1 axis ratios (Myers et al. 1996), 
which is geometrically closer to the sphere than the cylinder. So, the evolution lines of the 
two models in Figure |3] may have slightly faster infall speeds if the cores are spheroids, but 
the difference should be quite a bit less than 50% of the plotted value. 

Isothermality is another assumption in this study. Starless cores are observed to have 
temperature profiles decreasing inward because the interstellar radiation field heats the 
outskirts of the cores (e.g. Evans et al. 2001; Zucconi et al. 2001; Ward-Thompson et al. 
2002; Pagani et al. 2004; Shirley et al. 2005; Crapsi et al. 2007; Launhardt et al. 2013). We 
assumed isothermality because it simplifies the dynamics calculation by obviating the need 
to perform radiative transfer calculations. The effects due to non-isothermal temperature 
profile have been studied by Sipila, Harju, & Juvela (2011). They demonstrated that the 
critical sizes of BE-like cores depend on their temperature profiles; however, the critical 
size varies only about 6% (from 6.38 to 6.76) from the original critical size (6.5), while the 
temperature varies 35% (from 7 K to 9.5 K). Considering that the observed measurements 



-17- 

of ^max has a typical error of 50%, our results with an isothermal assumption are still 
applicable to real starless cores. 

Our analysis also ignores rotation of starless cores. Studies of molecular line widths 
toward starless cores show that cores slowly rotate (e.g. Heiles & Katz 1976; Arquilla & 
Goldsmith 1986; Goodman et al. 1993; Caselli et al. 2002). If the rotational energy is 
comparable to the gravitational energy, rotation cannot be ignored in studying dynamical 
stability of starless cores. The typical observed ratio of the rotational energy to the 
gravitational energy is about 0.02 (Goodman et al. 1993). So, our assumption of ignoring 
rotation is practical for typical starless cores. 

Magnetic fields may support starless core against gravitational collapse. Magnetic 
field strength of starless cores from Zeeman splitting measurement and from polarization 
observations range from 10 fiG to < 200 /iG (e.g. Ciolek & Basu 2000; Levin et al. 2001; 
Crutcher et al. 2004; Kirk, Ward- Thompson & Crutcher 2006; Turner & Heiles 2006). 
Magnetic pressure becomes comparable to the thermal pressure when the magnetic field 
strength reaches over 100 fiG. Our stability analysis is not accurate for starless cores 
with magnetic fields of > 100 fxG. Inclusion of magnetic term in the virial theorem is not 
practical for a dimensionless analysis because the inherent geometry of the field must be 
taken into account. Treatment of magnetic fields is beyond the scope of this current work. 

5.2. Comparison to Observed Starless Cores with Infall 

We compare the measured infall speeds and dimensionless sizes of observed starless 
cores which have clear infall signatures and well determined outer radii in Figure El The 
number of cores used in this study is limited because either only the density or the velocity 
structure is usually studied from observations of starless cores, but rarely both. The plotted 
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lines are the same with those in Figure [31 The observed starless cores are marked with 
black filled circles. The ^max value of L694-2 is from Harvey et al. (2003b). The peak 
infall velocity of L694-2 is measured by Lee et al. (2007) using the three hyperfine line 
of HCN J=l-0. The dimensionless size of L1544 is estimated by Kirk et al. (2005). The 
infall velocity of L1544 is studied by William, Lee & Myers (2006) with interferometric 
observation of N2H+ J=l-0. For L63, we quote the dimensionless size of ^max = 15 from 
Kirk et al. (2005). There is no published infall speed study of L63, therefore, we estimate 
the infall velocity by applying the two-layer model (equation (9) of Myers et al 1996) to 
the HCN lines observed by Sohn et al. (2007). From the HCN J=l-0 F=l-1 line, we 
estimate the infall speed of L63 to be about 0.1 km s^^. The dimensionless size of L1689B is 
measured by Dapp & Basu (2009) and the infall speed is measured by Bacmann & Pagani 
(2008). 

To normalize the infall speeds, internal temperature of those starless cores are required. 
Since internal temperatures of L63, L694-2 and L1689B are not studied with high precision, 
we deduced kinetic temperatures of the starless cores utilizing the hyperfine lines of NH3 
(1,1) and (2,2) inversion transitions. We fit ammonia spectra with five parameters including 
fisr, line width au, kinetic temperature Tk, optical depth of (1,1) transition ri, and filling 
fraction rjf using the same LTE assumption and method in Rosolowsky et al. (2008). 
We fit signals with Signal-to-Noise Ratio > 2 and estimate the goodness of fit using the 
reduced Xr- We calculated the reduced xl i^i the full five parameter space instead of using 
a covariant matrix to more accurately determine Icr uncertainties in the parameters. The 
results of the fittings are shown in Figure |H The black solid lines represent observed 
ammonia lines and red dashed lines represent fitted line profiles. Reduced Xr ^t the vicinity 
of the minimum value are plotted in Figure [51 The cross marks are the minimum points 
and the white lines represents Icr uncertainty space in T^ vs. Ti space at the minima of fisr, 
(7^, and rjf (see Table 2). The temperatures of L1544, L63, L694-2, L1689B from our study 
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are 8.84±0.30, 9.50±0.41, 9.14±0.2, and 12.0±1.7 K, respectively. Crapsi et al. (2007) 
also deduced average temperature of L1544 using the same transition of NH3 hyperfine 
lines, but with interferometric observations, to be 8.75K. Our results are consist within la 
uncertainty. The temperature of L694-2 also agrees with the previous finding of Williams 
et al. (2006) which is 9.3K. The temperature of L1689B has a relatively large uncertainty. 
This is because L1689B appears to be optically thin in the NH3 (1,1) line and determining 
its optical depth is degenerated with filling fraction. 

Combining all physical quantities, four starless cores are plotted in Figure O We mark 
error bars only for L694-2 because the quoted references do not provide any error for .^max 
and infall speeds for the other three cores. Generally, ^max has about 50% uncertainty, 
and errors of infall speeds are hard to estimate unless a starless core is mapped with high 
angular resolution and with multiple molecular lines. Accepting that the observed physical 
quantities for the starless cores quoted are the best determination to date, we may say 
that all four cores are located in the regime below the collapse of the critical BE sphere. 
Two cores, L1689B and L694-2, seem to collapse faster than the collapse of the critical 
uniform density sphere. These two starless cores may be undergoing enhanced collapse by 
external perturbations, for example, turbulence, or a sudden increase of external pressure. 
Observations of their surrounding environments are needed to understand the nature of the 
perturbation. The other two starless cores L1544 and L63 seem to be in quasi-equilibrium 
collapse. 

6. SUMMARY & CCONCLUSIONS 

Starless cores are the precursors of low-mass star formation. Their stability and 
dynamics are important issues in developing a comprehensive picture of how low-mass stars 
form. However, owing to the difficulties of observing internal motions of starless cores 
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through molecular line observations, only their density structures have been previously 
considered in investigating their stability and dynamics. In this study we find a stability 
condition for starless cores with internal motions. We first generalize the conventional 
stability condition of BE spheres by taking inward homologous motions into account. Even 
though homologous inward motions may not be a realistic internal motion that describes the 
observed internal motions of starless cores, it provides a rough estimate how the stability 
conditions of starless cores vary with the speeds of inward motions. Second, we suggest two 
limiting cases of spontaneous gravitational collapse: the collapse of the critical BE sphere 
and the collapse of uniform density spheres. The former model represents the slowest 
collapse by the self-gravity while the latter model sets the boundary for the fastest infall 
speed of spontaneous gravitational collapse at a given density. Any starless cores which 
considerably deviates from these two limiting models are likely to be strongly perturbed by 
the surrounding environments. 

We applied our study to four starless cores which have strong infall signatures (large 
blue peak asymmetry). Some cores appear to collapse faster than infall speeds permitted by 
spontaneous gravitational collapse. This result suggests that external perturbations may be 
important in the collapse dynamics of low-mass starless cores. Further study on the density 
and velocity structures of starless cores and their surrounding environments are required to 
obtain more general conclusions on the collapse dynamics of low-mass star formation. 

We thank anonymous referee for useful comments and corrections. 
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Table 1. Observed Starless Cores 



Source RA (J2000.0) DEC (J2000.0) I(T„ib) [K km/s] ai [Kkm/s] ffT„,b [K] 

L1544 05:04:17.2 +25:10:43.7 14.79 1.01 0.064 

L63 16:50:14.9 -18:06:22.5 10.22 0.72 0.107 

L694-2 19:41:04.3 +10:57:00.7 12.13 0.83 0.078 

L1689B 16:34:48.3 -24:38:03.5 2.962 0.26 0.093 



Table 2. Physical Quantities of Starless Cores Deduced from NH3 (1,1) and (2,2) 

Source v^g^, [kms^^] a^ [kms^^] Tk [k] ri rjf 

L1544 7.196 0.134+^05 8.84+°J° 3.17+°;^^ 0.80+^5 

L63 5.721 0.0825+nol 9.50+°;^} 2.62+0:1 0.80+°:°^ 

L694-2 9.595 0.127t°:|]°3 Q.IA+J'^f^ 2.361°;}^ 0.77t°:|]l 



L1689B 3.591 0.172^°;°^^ 12.01};^ 0.30+[]|[! 0.47+[j 



52 
0.21 



^Uncertainties of vis^ are not estimated 
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Uniform, P„,=cJ s({ J 
Uniform, P,„=0 
BE. P„,=cJsUJ 
BE. P„,=0 




Fig. 1. — The relationship between the inward homologous motion and the compression ratio 
for a uniform density core (black) and a Bonnor-Ebert core (blue). The pressure confined 
core (solid) is easier to compress than the core in a vacuum (dashed). Stronger compression 
occurs by a faster inward motion. 
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Fig. 2. — Critical size of the uniform density core (black) and the BE core (blue) with respect 
to the peak speed of inward motion. The solid and dashed lines denote cores bounded by 
two limits of the external pressure. 
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Fig. 3. — The evolution of the marginally unstable uniform density (red dashed) and BE 
(red solid) cores, and the critical sizes of the uniform density (black) and BE (blue) cores 
with respect to the peak speed of inward motion. The collapse of the BE core delivers the 
slowest infall motions, while the collapse of the uniform density core results in the fastest 
gravitational collapse. If a starless core goes through purely gravitational collapse, it is likely 
to be between the two evolution lines. If collapse of starless core is hindered by external 
perturbations, it may collapse slower than the collapse of the BE core. On the other hand, 
a starless core collapse faster than the collapse of the uniform density core is likely to be 
perturbed by external perturbations to enhance its infall speed. 
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Fig. 4. — Ammonia (1,1) and (2,2) lines observed toward starless cores. The black solid 
lines are the observed lines and the red dashed lines are fitted models. Fitting is done in 
the Gve parameters space including t>isr, line width, total optical depth of (1,1) transition, 
kinetic temperature, and filling factor. The best-fit parameters are listed in Table 2. 
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Fig. 5. — Reduced xl values in Tk vs. ti space. The other three parameters are fixed to he 
the best-fit values at the shown planes. The white crosses are the positions of the minimum 
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parameters are listed in Table 2. 
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Fig. 6. — Starless cores with clear infall signatures, L63, L1544, L1689B, and L694-2 are 
plotted with the generalized BE criterion and the limiting cases of the gravitational collapse 
for Figure 3. All cores are in unstable regime and gravitationally collapsing. The collapse of 
L694-2 and L1689B seem to be considerably enhanced by external perturbation. 



